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Abstract 

We reevaluate the two form factors relevant for the rfn second-class r decay mode, making systematic use of 
analyticity, unitarity, combined with updated inputs to the NLO chiral constraints. We focus in particular on the shape 
of the p resonance peak which is a background-free signature of a second-class current. Its dispersive construction 
requires the tjji — > nn scattering amplitude which we derive from a family of Khuri-Treiman equations solutions 
constrained with accurate recent results on the 77 — > 3;r Dalitz plot. 
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1. Introduction 

Weinberg remarked that exact isospin conservation 
by the strong interactions would imply selection rules 
for semi-leptonic weak decays [1]. Isospin conserva- 
tion, of course, is only approximate, being broken both 
in QCD because m„ nid and in QED because q,, + q^. 
Yet, no quantitative experimental evidence for these so- 
called second-class currents have been reported to date. 
Processes of this type in t decays are the rpi or ?/';r 
modes: because of parity conservation these decays 
must proceed through the 7=1 vector current which 
is even under a G-parity rotation, while r\n, rj'n are odd 
eigenstates of G-parity. The sensitivity of these modes 
to physics beyond the standard model (SM) has been 
discussed in refs \TJ3\. Within the standard model, mea- 
suring these T decay amplitudes would provide non triv- 
ial informations on matrix elements of the scalar opera- 
tor dfiuy^d and a related determination of the qq content 
of the scalar flo(980), flo(1450) resonances. The Babar 
collaboration has published upper bounds for the rjn, rj'n 
modes |i4j|5J. 

We reconsider here the theoretical expectation for the 
T — > riJiv mode in the SM. This problem was first ad- 



dressed in ref. [6|. We attempt to refine the theoret- 
ical prediction by exploiting the analyticity properties 
of the two form factors involved and combining them 
with chiral symmetry constraints. This proves particu- 
larly fruitful for the vector form factor, for which unitar- 
ity provides a simple relation with the isospin violating 
■q ^ 3>n decay amplitude. We will show how consid- 
erable recent progess in measuring this amplitude im- 
pacts the determination of the rjji vector form factor. It 
is usually expected that the integrated branching frac- 
tion should be dominated by the scalar rather than the 
vector form factor However, the main experimental ob- 
stacle to the observation of second-class amplitudes at 
B-factories is the pollution from first-class background 
contributions (e.g. from t"^ Tjji'^if'v, tjjt'^K^^ where 
the extra neutral particle escapes detection). While a 
peak in the rfii invariant mass at the fl()(980) mass may 
be present in background modes, a peak at the p(770) 
mass unmistakingly signals a second-class contribution. 

2. Dispersion relations and chiral symmetry 

The rjn form factors satisfy simple analyticity proper- 
ties, in exactly the same way as the more familiar nn or 
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Ktt form factors. They can be defined as analytic func- 
tions of the energy variable, s - (p^ + p„)^ with a right- 
hand cut and, furthermore, they are expected to behave 
as l/slog(i) when \s\ — > oo. Consequently, they satisfy 
unsubtracted dispersion relations (DR's). In practice, it 
is judicious to write DR's for the form factors multi- 
plied by weight functions like For instance, if we 
consider fTi^)! s^, the DR writes 



fTis) = fT(Q)+sfT{0)+ 



-r 



,disc[/fV)] 
as (1) 

Ami {s'Yis' - S) 



Thanks to the cutoff function, the integrand is dom- 
inated by the energy region below 1 GeV where we 
can evaluate the discontinuity, using unitarity, with only 
a few channels contributing (essentially, only a single 
channel). The price to pay is that we have to provide the 
values of the form factor and its derivative at s = 0. For 
this purpose, we can rely on three flavour chiral sym- 
metry, since the n and the 77 are both pseudo-Nambu- 
Goldstone bosons in this framework. 

The vector and scalar Tjn form factors are defined 
starting from the matrix element of the vector current 

(777r1«/£/|0)=- V2 [fTis) (p, - p„f+f""{s) (p, + p„r] 



and 



(2) 



At LO in the chiral expansion, the two form factors are 
constant and equal. 



V3(otj - OT„) 



: 0.99 X 10" 



There is no electromagnetic contribution at this order 
and the numerical estimate uses the chiral expansion of 
the mass differences - m^, - also at LO. 

The two form factors were computed at NLO in the 
chiral expansion by Neufeld and Rupertberger [7], in- 
cluding also the EM contributions at order e^. A re- 
markably simple expression emerges from their results, 
relating the value of the t/tt form factors at = to the 
ratio of the K^n'^ and K'^n^ form factors 



/f(0) 



1 

V5 



ft " (0) 



1 



log 



(3) 



Exploiting the recent results on K^^^ and TiTj^ decays from 
K factories (see e.g. f8l|) yields the most precise evalua- 
tion of the rjji form factors at s - 0, 



/f (0) = f^\0) = (1.49 ± 0.23) X 10-2 (4) 
which is significantly enhanced from its LO estimate. 



3. Vector form factor and i/ -» 3;r 

The discontinuity of fl" can be associated with a sum 
over intermediate states of the matrix element of the 
vector current 

lm{ijn^\uy'd\0) = ^ ^ t:^,^„MuJ^ d\0) . (5) 

n 

The derivation is valid in the unphysical situation where 
the T] meson is stable and we will assume than an ana- 
lytic continuation as a function of is possible. Be- 
low 1 GeV, the sum in eq. (j5]l is essentially saturated by 
the contribution of the lightest state n - n'^n^ (as it is 
strongly enhanced by its coupling to the p resonance). 
This leads to the following estimate of the discontinuity 
( for 4ml < ^ $ 1 GeV^) 



disc[/f (s)] = -0{s - Ami) x 
s - Ami 



(6) 



32n ^jA,j„(s) 



•I 



FUs) dzzT*{s,t{z)) 



with Ari„{s) - (s - ml){s - ml), m-t = m,, + m„. In this 
equation, Fy is the pion vector form factor, which is 
precisely known experimentally. In addition, one needs 
to evaluate the nn — > 777? amplitude projected on the P- 
wave, partly in an unphysical region (s < ml). It can be 
determined using its analyticity properties together with 
experimental constraints on ?;—> Stt decay. 

Combining analyticity with elastic unitarity for nn 
(re)scattering leads to a system of Khuri-Treiman (KT) 
equations ||9]. The solutions of these equations in their 
full generality were first discussed in refs. ifTOl [TTI . A 
subtle point, in particular, concerns the treatment of the 
singularities of the partial wave projected tttt — > t/jt am- 
plitude ( thus, T^^y^ ~ l/{s - ml)^l^ when s — > m^_) in 
the integrals. Eq. (j6| shows that these singularities af- 
fect also the computation of the vector form factor and 
must be treated by the same method. In practice, this 
leads to a distortion of the shape of the p resonance as 
compared to a naive vector meson dominance (VMD) 
approach. 

The authors of ref. IfTTII argue that a four-parameter 
family of solutions are relevant for the 77 — » 3?! de- 
cay, which we can also use for our problem. Assuming 
a sufficiently fast convergence of the three-flavour chi- 
ral expansion, they propose to determine all these four 
parameters by matching the dispersive and NLO chiral 
amplitude in a region around the Adler zero. Unfortu- 
nately, the amplitude obtained in this manner turns out 
not to be in agreement with the experimental results on 
the Dalitz plot parameters (see table[T]below). One must 
thus determine the KT solution parameters partly from 
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matching to the NLO amphtude and partly from fitting 
the experimental Dalitz plot data flTl . We perform here 
a fit analogous to ref. U2il but constraining the four KT 
parameters to be exactly real. From the NLO amplitude, 
we use the position of the Adler zero but not the value 
of the amplitude slope at this point. Thus, we do not 
attempt to determine the value of the quark mass ratio 
from the t] decay rate as in (TT\ (see also lfT3l [141) but 
take this ratio from the PDG (which leads to a value of 
the slope at the Adler zero differing from the NLO pre- 
diction by approximately 20%). 



param. 


experimental 


NLO Match. 


Fit 


a 
b 


-1.090 ± 0.005+[; ™^ 
0.124 + 0.006 + 0.010 


-1.300 
0.463 


-1.065 
0.159 


d 
f 


0.057 + 0.006![;™^ 
0.14 + 0.01+0.02 


0.069 
0.001 


0.066 
0.107 


a 


-0.0315 + 0.0015 


0.015 


-0.0355 



Table 1 : Comparison of the j] —t 37t Dalitz plot parameters obtained 
from KT solutions with experiment. The Dalitz parameters a, b, d, 
f refer to the charged decay mode and are taken from 1 15 |, while a 
refers to the neutral mode and the quoted value is taken from the PDG. 

The results of using a KT solution in the discontinuity 
relation (j6]l and then computing the form factor from the 
DR with weight functions l/s and l/s^ and using NLO 
chiral constraints like (|4]) is illustrated in fig.[T] The sta- 
bility with respect to the weight functions is satisfactory 
below 1 GeV. The results are also compared with the 
naive VMD where fl" would be simply proportional 
to the pion form factor normalized to the value Q at 
the origin. The dispersive calculation is seen to yield a 
significantly reduced resonance peak as compared to a 
naive VMD modelling. The influence of the KT param- 
eters is also rather significant and provide an idea of the 
uncertainties of this calculation. 



4. Scalar form factor model 

The scalar form factor (see Q) coincides with the Tjn 
matrix element of the derivative operator id^uy^^d and a 
discontinuity relation analogous to eq. (j5]l can be writ- 
ten. Since id^uyd is itself isospin suppressed as can be 
seen from the Ward identity idfjiy^d — {nid — mu)ud — 
eAfjUy'^dr, the sum over states can be resticted to 



In principle, one should include states \ny) in the sum as required 
by the second term in the Ward identity. These contributions as well 
as other EM ones where the photon line attaches to the charged lepton 
are neglected here. 



X 

O 



p' match one-sub 
p match two-sub 
improved one-sub 
improved two-sub 

/r(o) X ^'v'(^) 




0.8 1 1.2 

(GeV) 



Figure 1: Comparison of a naive VMD model for fl" (upper curve) 
with dispersive calculations using KT solutions. 



isospin conserving T-matrix elements T„^^„. Below 
one GeV, the state n - rjji dominates and we can write 

lmf^\s)^e(s-ml)x (7) 



16;rs 



For such a unitarity relation Watson's theorem implies 
that the phase of the form factor coincides with the rjn 
scattering phase shift in the region of elastic scattering. 
It is then natural to employ a phase dispersive represen- 
tation for the form factor, e.g. 



f^{s) = .C(0) 



(8) 



X exp 



( s(s - A,,^) 



ds'-r—r 



(m„+m^f s' {S - lS.,^„){s' - s) ^ 
A. 



(which uses the weight function lls{s - follow- 
ing lfT6l ). The values of the form factor at .? = and the 
Dashen-Weinstein point s - Ajj„ must be provided from 
NLO ChPT. A difficulty at this point is that the tjtt scat- 
tering phase shift is not measurable by the same meth- 
ods as used for nn or Kn. The experimental information 
concerns the properties of the resonances which couple 
to rjTT and the phase shift is constrained near the thresh- 
old by chiral symmetry. We will use a simple model 
proposed in ref. |17| which interpolates between these 
pieces of information. This model makes the plausi- 
ble prediction that the global features of r/n scattering 
are fairly similar to those of nn and nK scattering. The 
phase shift is a steadily raising function and inelasticity 
sets in rather sharply at a two-particle threshold (KK) 
close to a resonance (ao(980)). We can make use of 
this analogy to make a guess for the behaviour of the 
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Figure 2: Scalar form factor modulus from the phase representa- 
tion ^ illustrating the effect of varying sjjp . 

form factor phase in the inelastic energy region. In the 
case of nTT or nK the analogous phase can be determined 
by solving a system of coupled Muskhelishvili-Omnes 
equations using a known set of T-matrix elements. A 
common feature is that the phase drops sharply by ap- 
proximately 71 close to the inelastic threshold, which 
causes a dip in the modulus of the form factor. The exact 
point Sdip where this happens in the case of ijn cannot, 
of course, be known without actually solving the anal- 
ogous equations, but it seems plausible that this should 
be somewhere in between the two resonances ao(980) 
and flo(1450). If one of these resonances can be inter- 
preted as a tetraquark state (i.e. having a suppressed 
coupling to the M operator) then, from the analyticity 
point of vue, this corresponds to s^jp lying close to the 
corresponding resonance. Fig.|2]illustrates different val- 
ues for Sdip and fig.[3]shows the complete spectral func- 
tion for T — > rjTTv assuming Sdip to be midway between 
ao(980) and flo(1450). 

Finally, integrating over the spectral functions, we 
obtain the branching fractions shown in the last line 
of table |4] The results are preliminary and the errors 
quoted reflect only the uncertainty associated with the 
variation of Sdip- Our results tend to be in the lower 
range of previous evaluations. 
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Table 2: Our preliminary results for the central values of the Tin 
branching fraction: the vector and scalar contributions are shown sep- 
arately and compared to earlier evaluations. 
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